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ABSTRACT 

We are presenting a detailed parameter study of the time-dependent electron 
injection and kinematics and the self-consistent radiation transport in jets of in- 
termediate and low-frequency peaked BL Lac objects. Using a time-dependent, 
combined synchrotron-self- Compton and external- Compton jet model, we study 
the influence of variations of several essential model parameters, such as the elec- 
tron injection compactness, the relative contribution of synchrotron to external 
soft photons to the soft photon compactness, the electron-injection spectral in- 
dex, and the details of the time profiles of the electron injection episodes giving 
rise to flaring activity. In the analysis of our results, we focus on the expected 
X-ray spectral variability signatures in a region of parameter space particularly 
well suited to reproduce the broadband spectral energy distributions of inter- 
mediate and low-frequency peaked BL Lac objects. We demonstrate that SSC- 
and external-Compton dominated models for the 7-ray emission from blazars are 
producing significantly different signatures in the X-ray variability, in particular 
in the soft X-ray light curves and the spectral hysteresis at soft X-ray energies, 
which can be used as a powerful diagnostic to unveil the nature of the high-energy 
emission from BL Lac objects. 
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Subject headings: galaxies: active — galaxies: jets — radiation mechanisms: 
nonthermal 



1. Introduction 

The class of objects referred to as blazars consists of the most extreme examples of 
active galactic nuclei (AGNs), namely 7-ray loud, flat-spectrum radio quasars (FSRQs), and 
BL Lac objects. They have been observed in all wavelength bands — from radio through 
very-high energy (VHE) 7-ray frequencies. More than 65 blazars have been identified as 
sources of > 100 MeV emission detected by the EGRET telescope on board the Compton 
Gamma-Ray Observatory (CGRO) (Hartman et al. 1999), and at least 6 blazars have now 
been detected at VHE 7-rays (> 350 GeV) by ground-based air Cerenkov telescopes (for a 
recent review see, e.g., Buckley (2001)). Blazars exhibit variability at all wavelengths (von 
Montigny et al. 1995; Mukherjee et al. 1997) on time scales — in some cases — down to less 
than an hour (Gaidos et al. 1996). 

The broadband continuum spectra of blazars are dominated by non-thermal emission 
and consist of at least two clearly distinct, broad spectral components. A sequence of sub- 
classes of blazars can be defined through increasing peak frequencies and a decreasing dom- 
inance of the 7-ray output in terms of vF v peak flux along a sequence from FSRQs via 
low-frequency BL Lac objects (LBLs) to high-frequency peaked BL Lac objects (HBLs), 
which is also correlated with a decreasing inferred bolometric luminosity of the sources (Fos- 
sati et al. 1998). For recent reviews of the observational properties of blazars see, e.g., 
Sambruna (2000), Padovani & Urry (2001), or Bottcher (2002). 

Although all extragalactic sources detected by ground-based air Cerenkov telescope fa- 
cilities to date are HBLs, the steadily improving flux sensitivities and decreasing energy 
thresholds of those instruments provide a growing potential to extend their blazar source 
list towards intermediate and even low-frequency peaked BL Lac objects. The detection of 
such objects at energies ~ 40 - 100 GeV might provide an opportunity to probe the intrinsic 
high-energy cutoff of their spectral energy distributions (SEDs) since at those energies, 77 
absorption due to the intergalactic infrared background is expected to be negligible at red- 
shifts of z < 0.2 (de Jager & Stecker 2002). Such detections should significantly further our 
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understanding of the relevant radiation mechanisms responsible for the high-energy emission 
of blazars and the underlying particle acceleration mechanisms. 

The low-energy component of blazar SEDs is well understood as synchrotron emission 
from ultrarelativistic electrons in a relativistic jet directed at a small angle with respect 
to the line of sight. In the framework of leptonic models (for a review of the alternative 
class of hadronic jet models, see, e.g., Rachen (2000)), high-energy emission will result from 
Compton scattering of lower-frequency photons off the relativistic electrons. Possible target 
photon fields for Compton scattering are the synchrotron photons produced within the jet 
(the SSC process; Marscher & Gear (1985); Maraschi, Ghisellini, & Celotti (1992); Bloom & 
Marscher (1996)), or external photons (the EC process). Sources of external seed photons in- 
clude the UV - soft X-ray emission from the disk — either entering the jet directly (Dermer, 
Schlickeiser, & Mastichiadis 1992; Dermer & Schlickeiser 1993) or after reprocessing in the 
broad line region (BLR) or other circumnuclear material (Sikora, Begelman, & Rees 1994; 
Blandford & Levinson 1995; Dermer, Sturner, & Schlickeiser 1997) — , jet synchrotron radi- 
ation reflected at the BLR (Ghisellini & Madau 1996; Bednarek 1998; Bottcher & Dermer 
1998), or the infrared emission from circumnuclear dust (Blazejowski et al. 2000; Arbeiter, 
Pohl, & Schlickeiser 2002). 

According to the now well-established AGN unification scheme (Urry & Padovani 1995), 
blazars can be unified with other classes of AGN, in particular radio galaxies, through ori- 
entation effects. However, Sambruna, Maraschi, & Urry (1996) have pointed out that such 
orientation effects can not explain the differences between different blazar sub-classes. In- 
stead, it has been suggested that the sequence of spectral properties of blazars from HBLs 
via LBLs to FSRQs can be interpreted in terms of an increasing total power input into non- 
thermal electrons in the jet, accompanied by an increasing contribution of external photons 
to the seed photon field for Compton upscattering (Madejski 1998; Ghisellini et al. 1998). 
It has been suggested that this may be related to an evolutionary effect due to the gradual 
depletion of the circumnuclear material being accreted onto the central black hole (D'Elia 
& Cavaliere 2000; Cavaliere & D'Elia 2002; Bottcher & Dermer 2002). Detailed model- 
ing of blazars in the different sub-classes (FSRQs, LBLs and HBLs) seems to confirm this 
conjecture (for a recent review, see, e.g., Bottcher (2002)). 

As mentioned earlier, blazars tend to exhibit rapid flux and spectral variability. The 
variability is most dramatic and occurs on the shortest time scales at the high-energy ends 
of the two nonthermal spectral components of their broadband SEDs. Particularly interest- 
ing variability patterns could be observed at X-ray energies for those blazars whose X-ray 
emission is dominated by synchrotron emission. Observational studies of X-ray variability in 
blazars have so far focused on HBLs and, in particular, on the attempt to identify clear pat- 
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terns of time lags between hard and soft X-rays. However, such studies have yielded rather 
inconclusive and often contradictory results (e.g., for Mrk 421: Takahashi et al. (1996); Fos- 
sati et al. (2000); Takahashi et al. (2000); or PKS 2155-304: Chiapetti et al. (1999); Zheng 
et al. (1999); Kataoka et al. (2000); Edelson et al. (2001)). Instead, the so-called "spectral 
hysteresis" of blazar X-ray spectral variability may prove to be a more promising diagnostic 
of the physical nature of acceleration and cooling processes in blazar jets: When plotting 
the X-ray spectral hardness vs. the X-ray flux (hardness-intensity diagrams = HIDs), some 
HBLs (e.g., Mrk 421 and PKS 2155-304) have been observed to trace out characteristic, 
clockwise loops (Takahashi et al. (1996); Kataoka et al. (2000)). In terms of pure SSC jet 
models, such spectral hysteresis can be understood as the synchrotron radiation signature of 
gradual injection and/or acceleration of ultrarelativistic electrons into the emitting region, 
and subsequent radiative cooling (Kirk, Rieger, & Mastichiadis 1998; Georganopoulos & 
Marscher 1998; Kataoka et al. 2000; Kusunose, Takahara, & Li 2000; Li & Kusunose 2000). 
However, interestingly, such spectral hysteresis could not be confirmed in a recent series of 
XMM-Newton observations of Mrk 421 (Sembay et al. 2002). 

In LBLs, the soft X-ray emission is also sometimes dominated by the high-energy end of 
the synchrotron component (Tagliaferri et al. 2000; Ravasio et al. 2002), so similar spectral 
hysteresis phenomena should in principle be observable. However, those objects are gener- 
ally much fainter at X-ray energies than their high-frequency peaked counterparts, making 
the extraction of time-dependent spectral information an observationally very challenging 
task (see, e.g., Bottcher et al. (2002)), which may require the new generation of X-ray tele- 
scopes such as Chandra or XMM-Newton. Extracting the physical information contained in 
the rich X-ray variability patterns exhibited by BL Lac objects requires detailed theoreti- 
cal modeling of the time-dependent particle acceleration and radiation transport processes 
in the jets of blazars. Previous analyses of these processes (Kirk, Rieger, & Mastichiadis 
1998; Georganopoulos & Marscher 1998; Chiaberge & Ghisellini 1999; Kataoka et al. 2000; 
Kusunose, Takahara, & Li 2000; Li & Kusunose 2000; Krawczynski, Coppi, & Aharonian 
2002) have led to significant progress in our understanding of the particle acceleration and 
radiation mechanisms in HBLs, but were restricted to pure SSC models, with parameter 
choices specifically targeted towards HBLs. Consequently, those results may not be directly 
applicable to intermediate or low-frequency peaked BL Lac objects or even FSRQs. A no- 
table exception is a recent study by Sikora et al. (2001) (see also Moderski, Blazejowski, & 
Sikora (2002)), who included a significant contribution of external Compton radiation to the 
high-energy emission of blazars, and focused on the modeling of photon-energy dependent 
light curves and time lags between different frequency bands. They applied their results to 
the FSRQ 3C 279, and concluded that the correlated X-ray/7-ray variability of this quasar 
was inconsistent with X-rays and 7-rays being produced by the same radiation mechanism 
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because otherwise significant systematic time lags between the 7-ray and X-ray flaring be- 
haviour would be expected, contrary to the observations (e.g., Hartman et al. (2001)). 

In the present paper, we describe a newly developed combined SSC + ERC jet radiation 
transfer code, accounting for time-dependent particle acceleration and injection, radiative 
cooling, and escape, coupled to the self-consistent treatment of the relevant photon emission, 
absorption, and escape processes. In §2 we give a brief description of the underlying blazar jet 
model. The numerical procedure used in our code will be outlined in §3. We present results 
of a detailed parameter study, relevant for application to intermediate and low-frequency 
peaked BL Lac objects, in §4. We summarize in §5. 



2. Model Description 



The blazar model used for this study is a generic leptonic jet model. It is assumed 
that a population of ultrarelativistic, non-thermal electrons (and positrons) is injected at a 
generally time-dependent rate into a spherical emitting volume of co-moving radius Rb (the 
"blob"). The injected pair population is specified through an injection power L in j(t) and 
the spectral characteristics of the injected non-thermal electron distribution. We assume 
that electrons are injected with a single power-law distribution with low and high energy 
cutoffs 71 and 72, respectively, and a spectral index q so that the injection function Q e {^] t) 
[cm~ 3 s _1 ], in the co-moving frame of the emitting region, is 



Q i :Kr,t)=Q'o 3 (t) 1 ^ for 7 i<7<72 (1) 

with 

{ L inj ( f ) 2-q -r 1 q 

(2) 
V 6 'm e c2ln( 72 /7 1 ) 

where V b ' is the blob volume in the co-moving frame. 

The jet is powered by accretion of material onto a supermassive central object, which is 
accompanied by the formation of an accretion disk. For the purpose of this study, we have 
represented the disk by a standard Shakura-Sunyaev disk with a bolometric luminosity of 
L D = 10 45 ergs s _1 . The choice of this and several other standard parameters is motivated by 
a recent modeling study of the LBL W Comae (Bottcher, Mukherjee, & Reimer 2002). The 
randomly oriented magnetic field B is determined by an equipartition parameter e#, which 
is the fraction of the magnetic field energy density ub compared to its value for equipartition 
with the relativistic electron population in the emission region. As a consequence of this 
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parametrization, the magnetic field will gradually change throughout the evolution of the 
blob as particles are being injected and subsequently cool along the jet. The blob moves with 
relativistic speed v/c = fir = a/1 — 1/T 2 along the jet which is directed at an angle o b s (with 
/j, = cos# bs) with respect to the line of sight. The Doppler boosting of emission from the 
co-moving to the observer's frame is determined by the Doppler factor D — [F (1 — 

As the emission region moves outward along the jet, particles are continuously injected 
according to Eq. 1, are cooling, primarily due to radiative losses, and may leak out of the 
system. We parametrize particle escape through an energy-independent escape time scale 
^esc — V Rb/ C with r) > 1. This parametrization of particle escape can also be used to include 
adiabatic losses, which are not explicitly taken into account in our simulations. Radiation 
mechanisms included in our simulations are synchrotron emission, Compton upscattering of 
synchrotron photons (SSC = Synchrotron Self Compton scattering), and Compton upscat- 
tering of external photons (EC = External Compton scattering), including photons coming 
directly from the disk as well as re-processed photons from the broad line region. The broad 
line region is modelled as a spherical shell between rBLR,in = 0.2 pc and r BLR,out = 0.25 pc 
(we note, however, that the details of the radial distribution of the BLR material do not 
play an important role as long as rBLR.out — ^BLR,in ^$ ?"BLR,m), and a radial Thomson depth 
tt,blr which is considered a free parameter. 77 absorption and the corresponding pair pro- 
duction rates are taken into account self-consistently, using the general solution for the pair 
production rate of Bottcher & Schlickeiser (1997). However, in all simulations presented in 
this paper, the 77 opacity is <C 1 out to at least several tens of GeV so that 77 absorption 
and pair production do not play an important role. Motivated by the ~ 10 hr minimum 
variability time scale observed in W Comae (Tagliaferri et al. 2000), we choose R^ = 10 16 cm, 
and r = D — 10, which implies 8 ohs = 5.74°. 

Based on an equipartition parameter e# ~ 1, we expect typical magnetic field values 
of order B ~ 1 G (Bottcher, Mukherjee, & Reimer 2002), which implies a synchrotron 
cooling time scale (in the observer's frame) of electrons emitting synchrotron radiation at an 
observed energy E sy = 1 -Ekev keV of 

/ d \ - 3 / 2 / n \ - 1 / 2 
^°' 29 (tg) (io) ^h, (3) 

For X-ray photon energies, this is shorter than the dynamical time scale Rb/{Dc), in agree- 
ment with the approximately symmetric shape of the X-ray light curves generally observed 
both in W Comae (Tagliaferri et al. 2000) and BL Lacertae (Ravasio et al. 2002). 
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3. 



Numerical Procedure 



In order to treat the time-dependent electron dynamics and radiation transfer problem 
in the emitting volume, we solve simultaneously the kinetic equation for the relativistic 
electrons, 



Here, {dq / dt)\ oss is the radiative energy loss rate for the electrons, Q e (7,t) is the sum of the 
external injection rate Q 1 ^ from Eq. 1 and the intrinsic 77 pair production rate, %, hem (e, t) 
and n p h,abs(c, t) are the photon emission and absorption rates corresponding to the various 
radiation mechanisms, and t p h, C sc = (3/4) Rb/c. In Eq. 4, electron cooling is approximated 
as a continuous function of time (i.e., the energy of an individual electron is described as a 
differentiable function of time). This would be inaccurate if a significant contribution to the 
cooling rate were due to Compton scattering in the Klein-Nishina limit since in that case, the 
electron is transferring virtually all of its energy to a soft photon in a single scattering event. 
However, in the parameter ranges which we are primarily interested in, electron cooling is 
dominated by synchrotron losses and Compton scattering in the Thomson regime, for which 
Eq. 4 is a good approximation. In Eq. 5, the emissivity term contains the contribution 
from Compton scattering into a given photon energy interval. Since in all model situations 
considered here, the Thomson depth of the emitting region is tt *C 10~ 6 , the modification 
of the photon spectrum due to scattering of photons out of a given photon energy range is 
negligible. 

The relevant electron cooling rates and photon emissivities and opacities are evaluated 
using the well-tested subroutines of the jet radiation transfer code described in detail in 
Bottcher, Mause, & Schlickeiser (1997) and Bottcher & Bloom (2000). The full Klein-Nishina 
cross section for Compton scattering and the complete, analytical solution for the 77 pair 
production spectrum of Bottcher & Schlickeiser (1997) are used. The discretized electron 
continuity equation can be written in the form of a tri-diagonal matrix as in Chiaberge & 
Ghisellini (1999), which can be readily solved using the standard routine of Press et al. 
(1992). This procedure turns out to be very stable if, instead of the sharp cutoffs of the 
electron injection function (1), we introduce continuous transitions to very steep power-laws 
to mimic these cutoffs. Specifically, we add a low-energy branch with Q e {q;t) oc 7 2 for 





and for the photons, 



dn ph (e,t) 
dt 



^ph,cm 

(e,t) - n p h,abs(e,0 




(5) 
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7 < 7i> an( i Qe(l',t) oc 7 ^ +3 ^ for 7 > 72. After each electron time step, we update the 
photon distribution using a simple explicit forward integration of the discretized Eq. 5. 

We have carefully tested our code by running it with parameters identical to those used 
for Figs. 6 - 14 of Li & Kusunose (2000), who are using a very similar numerical approach. 
We find very good agreement with their results, with only minor discrepancies which are 
due to our replacing the high-and low-energy cutoff of the electron injection spectrum by 
continuous transitions to very steep power-laws as described above. Specifically, this results 
in more moderate spectral indices of the synchrotron spectra just beyond the high-energy 
cutoffs. 



4. Numerical Results 

We have performed a large number of simulations, studying the influence of various 
model parameters on the resulting broadband spectra, light curves, and X-ray hardness- 
intensity diagram tracks. In each one of our simulations, we have assumed an underlying, 
quiescent injection power of L?^ = 10 38 ergs s _1 , on top of which we inject particles with 
various flaring injection powers in the range 10 40 ergs s" 1 < Lf n j < 10 43 ergs s -1 . As a 
standard model setup, we choose an injection electron function given by 71 = 10 3 , 72 = 10 5 , 
and q = 2.5. In the base model, we have Lf n - = 10 41 ergs s -1 , extending as a step function 
in time over 2 dynamical time scales, At[ n j = 2R b /c in the co-moving frame. The injection 
event is centered around a distance xq = 0.1 pc from the central engine. The BLR Thomson 
depth is chosen to be in the base model. 

Subsequently, we investigate the influence of changing (1) the flaring injection power, 
(2) the BLR Thomson depth and, accordingly, the contribution of external photons to the 
soft photon field for Compton scattering, (3) the electron injection spectral index q, (4) the 
duration of the flaring injection event, (5) the time profile of the flaring injection event, (6) 
the electron escape time scale parameter. The parameters used for the individual runs are 
quoted in Table 1. 

Our base model is simulation no. 2. In Figs. 1 and 2, we have compiled a sequence 
of co-moving electron spectra, snap-shot SEDs, and the time averaged photon spectrum 
resulting from our base model. Light curves at 3 selected X-ray energies as well as in the 
optical (R-band) and at hard X-rays (30 keV) are plotted in Fig. 3, and tracks in the 
hardness-intensity diagrams (HIDs) at three different X-ray energies are compiled in Fig. 
4. The figures illustrate the gradual build-up of the electron density in the emission region, 
competing with radiative cooling, which is faster than the injection time scale at electron 
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energies of 7 > 10 4 . Radiative cooling is the dominant process affecting the electron spectra 
after the end of the flaring injection episode at t — 2 t dyn . The time-dependent photon spectra 
as well as the light curves demonstrate that we do not expect significant peak time delays 
within the synchrotron component at frequencies v > 10 14 Hz, but that the high-energy 
(SSC) component is delayed by ~ 1 dynamical time scale due to the gradual accumulation 
of seed photons for Compton scattering. The figure also indicates the very moderate flux 
variability at energies just above the synchrotron cut-off, which is located at ~ 1 keV in 
our example. Fig. 4 illustrates the spectral hysteresis phenomenon (keeping in mind that 
additional contributions from previous injection episodes should close the tracks in the sense 
that they are expected to start out near the end points of the tracks shown in the figure). 
In agreement with Li & Kusunose (2000) we find that — at least in this generic case — the 
spectral hysteresis tracks can change their orientation from clockwise to counterclockwise 
as one goes from photon energies below the synchrotron cutoff to energies above the cutoff, 
where the spectrum is dominated by Compton scattering (SSC). 

In the following, we are focusing on the time-averaged photon spectra, the light curves, 
and the X-ray spectral hysteresis, and investigate how those aspects are affected by variations 
of individual parameters. 

4.1. Electron Injection Power 

The effect of an increasing injection power - - corresponding to a higher density of 
injected, relativistic particles in the emitting region — is illustrated in Figs. 5-7. In 
addition to a corresponding increase in the overall bolometric luminosity, this leads also 
to a stronger relative energy output in the SSC-dominated Compton emission at X- and 
7-ray energies, as expressed, e.g., in Eq. (19) of Chiang & Bottcher (2002). The photon 
spectral index of the time-averaged emission at optical - soft X-ray frequencies remains 
robust at a Q _ x ~ 1.25 due to optically thin synchrotron emission from the cooled electron 
spectrum with injection spectral index q = 2.5. As the bolometric luminosity (and the 
electron cooling) becomes dominated by the SSC mechanism, one would expect that this 
changes to the canonical value of a D -x = 1-5, which is a result of the decaying electron 
cooling rate in an SSC-dominated cooling scenario (Chiang & Bottcher 2002). However, in 
the situation simulated here, the first-order SSC peak is rapidly (within ~ 3 tdyn) decaying 
into the keV energy range and dominating over the instantaneous synchrotron emission at 
UV - X-ray energies. This leads to a significant hardening of the time-averaged optical - 
X-ray spectrum, which even becomes inverted in vF v space in our most extreme test case 
(simulation no. 8). 
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In the light curves (see Fig. 6), the more rapid electron cooling with increasing electron 
injection power manifests itself in an overall increasing amplitude of variability at all energies. 
In particular, as SSC cooling becomes more important, even the harder X-rays begin to 
exhibit significant variability on the dynamical time scale, in contrast to the synchrotron- 
cooling dominated cases. Furthermore, while for very low injection powers, the synchrotron 
cooling time scale for optical synchrotron emission is comparable to the injection time scale, 
resulting in a time delay of ~ a few hr between X-ray and optical emission, the optical light 
curve peaks at the end of the injection episode for higher injection powers, simultaneously 
with the X-rays. However, when SSC cooling becomes dominant, the gradually increasing 
energy density in the soft photon field during the injection episode actually has the effect 
that the X-ray light curves are peaking at the beginning of the injection episode, which 
would, again, lead to a time delay of ~ a few hr between X-ray and optical flares. 

Fig. 7 illustrates how the tracks in the HIDs at different X-ray energies are drastically 
changing for different injection powers. In particular at X-ray energies just below or at the 
synchrotron cutoff (~ 1 keV), the flux maxima are occurring at significantly different values 
of the local spectral index a for different values of the injection power. Specifically, the 
local spectral indices at the time of the peak flux are significantly smaller (harder) for larger 
values of the injection power. Obvious changes in the orientation of the spectral hysteresis 
tracks are not found in these simulations. 

As mentioned earlier (see Eq. 3), in the test cases investigated here, the synchrotron 
cooling time scale of electrons emitting synchrotron radiation at X-ray energies, is shorter 
than the dynamical time scale, which is of the same order as the injection time scale. Conse- 
quently, our results can be qualitatively compared to those of Li & Kusunose (2000) for the 
"short cooling time limit" , bearing in mind that the parameter values in our simulations have 
been chosen appropriate for intermediate and LBLs, while Li & Kusunose (2000) focused on 
the application to the HBL Mrk 421. In particular the light curves displayed in their Figs. 
9-11 exhibit the same general trends as we have found in our set of simulations. 

4.2. External Photons 

In order to investigate the influence of an increasing contribution of external photons 
to the soft photon field for Compton scattering, we performed a series of simulations with 
increasing values of Tt,blr, from to 1. We note that the contribution from direct accretion 
disk photons to the photon energy density in the emitting region is negligible in our base 
model. Thus, in the following, the external photons are primarily accretion disk photons 
reprocessed in the BLR. An additional component due to direct accretion disk photons 
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would become energetically important only for injection flares significantly contributing at 
< 10- 2 L^ 2 /[(S/G)r!] pc, where L 45 is the accretion disk luminosity in units of 
10 45 ergs s _1 and I\ = T/10. For a recent discussion of this radiation field on the high- 
energy light curves see, e.g., Dermer & Schlickeiser (2002). 

The time-averaged photon spectra from our simulations with varying BLR Thomson 
depth are shown in Fig. 8, which clearly shows the emergence of the external Compton 
(EC) component at GeV 7-ray energies. The impact of this additional emission component 
on the lower-frequency emission is small as long as its bolometric energy output is smaller 
than or comparable to the synchrotron vF v flux. Only when EC cooling becomes dominant 
over synchrotron cooling, are the effects on the synchrotron + SSC dominated portion of the 
spectrum (radio frequencies - MeV 7-rays) noticeable. This is the case when the comoving 
energy densities of the magnetic field and the external photons, u' B and u' ext , respectively, 
become comparable. From 



2 



L45 Tt,BLR I\ 



- where ro.2 is the radius of the inner boundary of the BLR in units of 0.2 pc — we can esti- 
mate that this happens at tt,blr ~ 0.1. Specifically, for tt,blr ^ 0.1, the effect of dominant 
external-Compton cooling results in a reduction of the time-averaged emission around the 
synchrotron peak — leading to a spectral hardening of the synchrotron emission — , and a 
suppression of the SSC emission. The suppression of the low-frequency synchrotron emission 
can be explained as the combined effect of two causes: First, electrons at energies above the 
low-energy cut-off are radiatively cooling on a time scale much shorter than the dynamical 
one (see Eq. 7 below). Consequently, the particle spectrum of high-energy electrons injected 
during the flare is rapidly depleted within less than one dynamical time scale from the end 
of the flaring episode. However, after this episode, we are still injecting electrons (although 
at a much smaller rate corresponding to L^) into the blob, which continue to present a 
high-energy electron population. In the case of the extremely fast cooling rate at tt,blr ~ 1, 
this additional high-energy electron population leads to a significant additional contribution 
of synchrotron emission at intermediate energies beyond the down-shifted maximum energy 
of electrons injected during the flare and, consequently, to a hardening of the synchrotron 
spectrum. 

The second cause of spectral hardening of the low-frequency synchrotron spectrum is 
related to our parametrization of the magnetic field in terms of an equipartition parameter 
€3- For electron injection spectral indices q > 2, most of the energy in the electron population 
is carried by electrons near the low-energy cutoff. Consequently, the co-moving magnetic 
field will decay on a time scale given by the radiative cooling time scale of the lowest-energy 
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electrons, which is 

t' B = *ec(7i) « 4 x 10 4 r } 2 s (7) 

r T,BLR7l,3 -^45 1 1 

where 7^3 = 71/IO 3 . For 7t,blr ^ 0.1, the time scale (7) is of the order of or shorter than the 
dynamical time scale. In contrast, the energy density of the external radiation field remains 
approximately constant during the evolution of the blob. Consequently, as radiative cooling 
proceeds and shifts the electron distribution towards lower energies, a steadily decreasing 
fraction of the electron energy will be converted to synchrotron radiation. This leads to a 
spectral hardening of the time-averaged synchrotron spectrum compared to the canonical 
a = q/2 spectrum above the break frequency v\ = VLfilj °f a cooling electron population 
in a constant magnetic field B = B with vlq = eB / '(2nm e c) . At frequencies below u\, 
the time-averaged spectrum would have a slope of a = 1/2 in the constant-magnetic- field 
case. The steepening of this low-frequency part of the synchrotron spectrum in the case of 
a magnetic field proportional to the equipartition value can be derived analytically in the 
following way. 

At late times, flaring electron distribution has basically collapsed to a 8 function in 
electron energy, i. e. n(7 e ,t) oc 5{j e — 7[t]). Then, the synchrotron emission coefficient will 
behave as 



ju tB y(t) OC % y 5(u -7[t] 2 Z/ L [t]) 

oc >y[t] 3 6(v->y[t] 5 ' 2 v L ,o). 



(8) 
(9) 



where 7 sy oc 7(t) 3 is the synchrotron loss rate. The additional factor of 7 results from the 
equipartition prescription, B 2 /8n = eBn e ^f. The electron evolution is governed by EC losses, 
and since the EC photon field is essentially constant over the flare episode, the electron 
Lorentz factor still evolves according to 7 — (A/3)(a T /m e c)u' cxt r ) 2 . The time-integrated 
synchrotron spectrum is then 



(jv,syn)t 



OC 



OC 



J ^7M 3 <^-7M 5/ V,o) 



l[t\ 



dt 



OC V 



(10) 
(11) 



where we have used v oc 7 5 / 2 . The short-dashed curve in Fig. 8 shows the result of a test 
calculation in which we held the magnetic field constant, while all other parameters were 
identical to the tt,blr — 1 simulation, in order to verify that the spectral hardening at radio 
frequencies is indeed partially a consequence of our magnetic field parametrization. 
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Fig. 9 shows that the impact of a strong external Compton component on the optical 
and X-ray light curves is very moderate. In particular, the light curves at energies below 
the synchrotron cutoff remain virtually unchanged. However, a strong external Compton 
component, dominating the bolo metric luminosity (t^blr ^ 0.1 in our case) leads to a 
significantly faster decay of the light curves at X-ray energies beyond the synchrotron cutoff. 

Just as the light curves, also the X-ray spectral hysteresis characteristics remain virtually 
unchanged, even in the case of a strongly dominant external Compton component, except 
for a moderate softening of the local spectrum in the decaying phase of the flare at energies 
beyond the synchrotron cutoff (see Fig. 10). 

4.3. Electron Spectral Index 

The value of the electron injection spectral index q should be rather easily determined by 
measuring the time-averaged optical - X-ray spectral index of the strongly cooled synchrotron 
spectrum. As illustrated in Fig. 11, this spectral change is accompanied by a shift of the 
SSC peak towards higher frequencies as the injection spectrum hardens: For q > 2.5, the 
SSC peak is located at esse ~ 7i e s y , while for q < 2.5, it shifts towards esse ~ l2 e sy As 
illustrated in Fig. 12, the characteristics of the light curves are only marginally different for 
different values of q. The X-ray spectral hysteresis at X-ray energies below the synchrotron 
cutoff is obviously shifted according to the change in injection spectral index, but its basic 
characteristics remain unchanged (see Fig. 13). An interesting qualitative change of the 
spectral hysteresis can be seen for energies just above the synchrotron cutoff: While for hard 
electron injection spectra (q < 2.5), the peak flux is reached at a steep local X-ray spectrum 
(i.e. the spectrum is dominated by the synchrotron component at the time of peak flux), a 
soft injection spectrum (g > 2.5) leads to a hard local spectral index at the time of peak flux 
(i.e. the spectrum is dominated by the SSC component at that time). 

4.4. Duration of the Flare 

In order to investigate the influence of the duration of the electron injection event 
causing the flare, we have performed simulations with 3 different values of t m j, where we 
kept the total energy input during the flare constant (simulations no. 2, 11, and 12). The 
time-averaged SEDs from those simulations are virtually indistinguishable. Fig. 14 shows 
that the longer injection time scale at a lower injection power leads to a more gradual rise 
of the light curves at energies above the synchrotron cut-off, whereas at lower frequencies, 
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this leads to a rather marginal modification of the rising portion of the light curve, followed 
by an extended plateau until the end of the injection episode. The decaying portion of 
the light curves at optical and X-ray frequencies seems to be virtually independent of the 
duration of the injection event. In the X-ray spectral hysteresis (Fig. 15), we find that 
for longer electron injection events, the rising-flux portion of the track in the HID at X-ray 
energies below the synchrotron cutoff occurs with increasingly softer local spectra, whereas 
the decaying portion remains virtually unchanged. At energies above the synchrotron peak, 
the trend concerning the rising portion is opposite: With longer duration of the injection 
episode, the local spectra are becoming harder. 

From an observational point of view, one should be able to distinguish situations in 
which the duration of the flaring event is substantially longer than the dynamical time scale, 
by virtue of the extended high-flux plateaus at energies below the synchrotron peak, which 
are generally not observed in BL Lac objects in which there is evidence for a substantial 
contribution from synchrotron emission to the X-ray flux (e.g., Tagliaferri et al. (2000); 
Ravasio et al. (2002)). This seems to indicate that the generic situation of ii n j ~ 2td yn might 
be a realistic assumption. 



4.5. Time Profile of the Flare 

The step function time profile of the electron injection power during the flare is certainly 
a rather crude over-simplification of any realistic acceleration scenario. In order to investigate 
whether this particular choice of the time profile has a significant impact on our results, we 
have calculated an additional set of simulations (nos. 13 - 15), with triangular injection 
profiles. Here, we have introduced a linear rise and decay of the injection power on time 
scales t r and td, respectively, and have chosen the maximum of the profile to be twice the 
injection power of the step-function case in order to keep the total injected energy at the 
same value as in our base model. 

The time-averaged photon spectra for all of these cases are virtually identical. Equally, 
the light curves at energies above the synchrotron cutoff are only marginally affected by the 
details of the injection time profile (see Fig. 16), while at lower energies the light curves tend 
to track the injection time profile to a certain extent during the rising portion of the light 
curves. The decaying portion of the light curves is generally independent of the injection time 
profile for all optical and X-ray photon energies. Fig. 17 illustrates that the impact of the 
detailed injection profile on the X-ray spectral hysteresis characteristics is rather moderate, 
even for the extreme (and very artificial) cases of the time profiles c (gradual rise over 2tdyn, 
sharp decay) and d (sharp rise, gradual decay over 2t dyn ) illustrated here. Its influence is 
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basically restricted to the flux-rise portion of the HID track and to photon energies below the 
synchrotron cut-off, where the local spectra tend to be harder for time profiles with maxima 
closer to the onset of the flare. The HID tracks for all our test cases show almost identical 
spectral indices at the time of maximum flux. 



At the low-energy end of the electron spectra, particles cooling down from higher energies 
will either accumulate and build up a 7~ 2 power-law spectrum at energies below 71, or escape, 
depending on the value of the escape time scale parameter 77. We can define a critical escape 
parameter 



for which the escape time scale for electrons at energy 72 equals the synchrotron cooling time 
scale. For rj < i] cr , the electron distribution will maintain a sharp low-energy cutoff at 71, 
while for 77 ^> i] CT , a -2 low-energy power law will develop. In order to investigate whether 
our choice of r\ = 10 has a significant impact on our results, we have done test simulations 
with 77 = 3 and r\ = 30, respectively. We find all relevant aspects — the time-averaged 
spectra, the monoenergetic light curves, and the spectral hysteresis curves — to be virtually 
independent of r\ within reasonable bounds. 



In the previous subsections, we have discussed the impact of various parameters on the 
broadband SED and the X-ray variability characteristics of models for intermediate and low- 
frequency peaked BL Lac objects. There are still a few more parameters left which we have 
fixed in our suite of test simulations. In particular, the choice of the cutoffs of the electron 
distribution, 7! and 72, the magnetic- field equipartition parameter e#, the Doppler factor 
D, and the size of the emitting region, Rb, may have an impact on the spectral variability 
characteristics. However, those parameters can generally be constrained rather well through 
the observed overall spectral characteristics, and through variability time scale considerations 
(see, e.g., Tavecchio, Maraschi, & Ghisellini (1998); Costamante & Ghisellini (2002)). The 
values adopted here are representative for typical LBLs like BL Lacertae (Madejski et al. 
1999; Bottcher & Bloom 2000) or W Comae (Tagliaferri et al. 2000; Bottcher, Mukherjee, 
& Reimer 2002). Furthermore, significantly different values of 72, e B) and D would shift the 



4.6. Electron Escape Time Scale 




(12) 



4.7. Other parameters 
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synchrotron cutoff out of the X-ray regime so that the diagnostics developed here may not 
be applicable to the X-ray variability of BL Lac objects, which is the focus of this paper. 

5. Summary and Conclusions 

We have presented a detailed parameter study of the time-dependent electron injection 
and kinematics and the self-consistent radiation transport in jets of intermediate and low- 
frequency peaked BL Lac objects. Those objects are currently of great interest as the 
steadily improving capabilities of current and future air Cerenkov detector facilities might 
allow the detection of this class of blazars at multi-GeV energies in the near future. At 
the same time, some of these objects exhibit interesting X-ray variability features which can 
now be studied in detail with the new generation of X-ray telescopes, in particular Chandra 
and XMM-Newton. Furthermore, the GLAST mission, scheduled for launch in 2006, is 
expected to detect many more BL Lac objects at multi-MeV - GeV energies and bridge the 
observational gap between the energy ranges previously covered by the EGRET instrument 
on board CGRO, and the ground-based air Cerenkov facilities. 

In our study, we have focused on the impact of various specific parameter choices and 
variations on the broadband SEDs, optical and X-ray light curves, and the spectral hysteresis 
phenomena previously observed in several high-frequency peaked BL Lac objects, but also 
expected to be observable in intermediate and low-frequency peaked BL Lacs. 

Very important conclusions can be drawn from a comparison of our results concerning 7- 
ray bright sources dominated by either SSC or external Compton emission. For a given level 
of flux at GeV energies, at a level comprable to or exceeding the synchrotron vF v peak flux, 
those two scenarios should be clearly distinguishable by virtue of the X-ray variability during 
flaring episodes: If roughly symmetric flare time profiles at soft X-rays below the synchrotron 
cutoff are observed, and the local, time-resolved X-ray spectra are soft — consistent with a 
case with negligible 7-ray emission — the 7-ray emission might be dominated by external- 
Compton emission (see Figs. 8 - 10). However, if the X-ray time profiles show a clear sign of a 
very rapid rise and more gradual decay, and the time-resolved spectra are significantly harder 
than corresponding to a case without strong 7-ray emission, we expect a strong contribution 
from the SSC mechanism to the 7-ray emission (see Figs. 5-7). Most notably, this diagnostic 
of the 7-ray emission mechanism does not require a detailed spectral measurement — which 
will be hard to achieve and require long integration times, even with GLAST — , but only 
a rough estimate of the GeV flux. It relies primarily on X-ray variability studies, combined 
with the type of detailed time-dependent radiation modeling presented in this paper. 



-17- 



We have shown that our results are not significantly impacted by the special choice 
of poorly determined parameters like the details (exact duration and time profile) of the 
acceleration / injection events leading to flaring activity, or the electron escape time scale 
parameter. Consequently, the X-ray variability of the high-frequency end of the synchrotron 
emission in intermediate and low-frequency peaked BL Lac objects can be used as a very 
robust diagnostic to unveil the nature of the high-energy emission in this type of blazars. 

We thank C. D. Dermer for discussions and very detailed and useful comments, and 
the anonymous referee for a careful review and a helpful report. The work of MB was 
supported by NASA through Chandra Postdoctoral Fellowship Award no. 9-10007, issued by 
the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory 
for and on behalf of NASA under contract NAS 8-39073. 
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Fig. 1. — Sequence of co-moving electron spectra in the emission region for our base model, 
simulation no. 2. The curves are labeled by time in multiples of the dynamical time scale, 
tjs = R b /(D C ) = 3.33 x 10 4 s. For model parameters, see Tab. 1. 
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v[Hz] 

Fig. 2. — Sequence of spectral energy distributions (upper panel) and the time averaged SED 
(lower panel) from our base model, simulation no. 2. The dotted vertical lines indicate the 
frequencies at which light curves and hardness-intensity correlations have been extracted; 
the dashed vertical lines indicate the remaining two frequencies at which light curves have 
been extracted. 
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Fig. 3. — Light curves at three different X-ray energies, in the optical (R-band) and at 
hard X-rays for our base model. Upper panel: Absolute flux values; lower panel: fluxes 
normalized to their peak values. The vertical dotted line indicates the end of the flaring 
injection episode. 



-24- 




Fig. 4. — Tracks of the simulated spectra in the hardness-intensity diagrams at three different 
X-ray energies for our base model, a is the local energy spectral index F v oc v~ a at the 
respective photon energy. Stars indicate the locations at multiples of the dynamical time 
scale during the simulation. 
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Fig. 5. — Time averaged photon spectra for different values of the flaring injection power Lf n ^, 
from simulations no. 1, 2, 3, and 8. The dotted vertical lines indicate the frequencies at which 
light curves and hardness-intensity correlations have been extracted; the dashed vertical lines 
indicate the remaining two frequencies at which light curves have been extracted. 
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Fig. 6. — Optical and X-ray light curves for different values of the flaring injection power 
Lf-, from simulations no. 1, 2, and 3. 



-27- 



0.7 |! 0.3 keV 



-0.3 




F v /F v 



max 



Fv/Fv 



max 



Fv/Fv 



max 



Fig. 7. — Tracks in the harness-intensity diagrams at X-ray energies for different values of 
the flaring injection power Lf n -, from simulations no. 1, 2, 3, and 8. Symbols are the same as 
in Fig. 5: short-dashed (10 40 erg/s), solid (10 41 erg/s), dot-dashed (10 42 erg/s), long-dashed 
(10 43 erg/s). Stars indicate the locations at multiples of the dynamical time scale during the 
respective simulations. 
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Fig. 8. — Time averaged photon spectra for different intensities of the external soft radiation 
field, parametrized through different values of t T blR) from simulations no. 2, 5, 6, 7, and 18. 
The short-dashed curve shows a test simulation with constant magnetic field, while all other 
parameters were identical to the tt,blr = 1 simulation, illustrating the effect of our magnetic 
field parametrization as a constant fraction of the equipartition magnetic field. The dotted 
vertical lines indicate the frequencies at which light curves and hardness-intensity correlations 
have been extracted; the dashed vertical lines indicate the remaining two frequencies at which 
light curves have been extracted. 
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Fig. 9. — Optical and X-ray light curves for different intensities of the external soft radiation 
field, parametrized through different values of tt,blr, from simulations no. 2, 6, and 7. 
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Fig. 10. — Tracks in the harness-intensity diagrams at X-ray energies for different intensities 
of the external soft radiation field, parametrized through different values of tt,blr, from 
simulations no. 2, 6, and 7: solid (tt,blr = 0), dashed (tt,blr = 0.1), dot-dashed (tt,blr = 
1). Stars indicate the locations at multiples of the dynamical time scale during the respective 
simulations. 
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Fig. 11. — Time averaged photon spectra for different values of the electron injection spec- 
tral index q, from simulations no. 2, 9, and 10. The dotted vertical lines indicate the 
frequencies at which light curves and hardness-intensity correlations have been extracted; 
the dashed vertical lines indicate the remaining two frequencies at which light curves have 
been extracted. 



- 32 - 




R-band 

0.3 keV 

--- 1 keV 
- - 5keV 
30 keV 



1 2 3 4 5 6 
t ob s[10 5 s] 

Fig. 12. — Optical and X-ray light curves for different values of the electron injection spectral 
indes q, from simulations no. 2, 9, and 10. 
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Fig. 13. — Tracks in the harness-intensity diagrams at X-ray energies for different values of 
the electron injection spectral index q, from simulations no. 2, 9, and 10: dashed (q = 2.2), 
solid (q = 2.5), dot-dashed (q = 2.8). Stars indicate the locations at multiples of the 
dynamical time scale during the respective simulations. 
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Fig. 14. — Optical and X-ray light curves for different values of the duration of the flar- 
ing injection time scale, keeping the total energy input during the injection event constant 
(simulations no. 2, 11, and 12). 
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Fig. 15. — Tracks in the harness-intensity diagrams at X-ray energies for different values 
of the flaring injection time scale, keeping the total energy input during the injection event 
constant (simulations no. 2, 11, and 12): solid (t^ = 2td yn ), dashed (ti n j = 4td yn ), dot- 
dashed (t in j = 6tdyn)- Stars indicate the locations at multiples of the dynamical time scale 
during the respective simulations. 
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Fig. 16. — Optical and X-ray light curves for different time profiles of the electron injection 
power, keeping the total energy input during the injection event constant (simulations no. 
2, 14, and 15). The small insets illustrate the injection time profiles (L in j vs. time). 
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Fig. 17. — Tracks in the harness-intensity diagrams at X-ray energies for different time 
profiles of the electron injection power, keeping the total energy input during the injection 
event constant (simulations no. 2, 14, and 15): solid (profile a = step function), dashed 
(profile c = triangular profile with t r = 2td yn and td — 0), dot-dashed (profile d = triangular 
profile with t r — and t d = 2t dyn ). Stars indicate the locations at multiples of the dynamical 
time scale during the respective simulations. 
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Table 1. Parameters of the simulations used for our parameter study. Throughout the 
series of simulations, we have used L?^ = 10 38 ergs s -1 (electron injection luminosity during 
quiescence), 71 = 10 3 (low-energy cutoff of injected electron spectrum), 72 = 10 5 
(high-energy cutoff of injected electron spectrum), D = 10 (Doppler boosting factor), 
Rb = 10 16 cm (blob radius), and e# — 1 (magnetic-field equipartition parameter). The time 
profiles in the 6th column are: (a) step function, (b) triangular with linear rise and decay 
with equal time scales (t r ^ = U n j/2), (c) linear rise and instantaneous drop (t T = £ m j), (d) 
instantaneous rise and linear decay (t d = t in j). For time profiles (b) - (d), the parameter 
L^j is the maximum injection power, r] is the electron escape time parameter, defined by 
£ e ,esc = ^Rb/c. The parameters changed with respect to the base model (no. 2) printed in 

boldface. 
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2.5 


2 idyn 


a 


10 


2 


10 41 
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*In this simulation, the magnetic field was held constant at B = 0.4 G. 



